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1. Introduction 

This paper is motivated by questions about the behavior of solutions of elliptic and para- 
bolic equations with low regularity drift terms. A classical example is 

d t u + b ■ Vu - Au = (1.1) 

considered in K n x]0, oo[, where b is a time-dependent vector field in MJ 1 . Of particular interest 
to us will be the case of divergence-free b (i.e., divfr = 0), which is relevant for applications 
to incompressible flows. 

To describe the regularity conditions on the drift term, it is useful to recall some elementary 
dimensional analysis. Equation fll.ip is invariant under the following scaling transformations: 

u(x, t) ->■ u w (x, t) = u(Xx, \H) , (1.2) 
b(x, t) ->■ b w (x, t) = \b(\x, \H) , (1.3) 

where A > 0. Following the usual convention (see, e.g., |3j), we can say that u has dimension 
and b has dimension — 1. The classical theory (see, e.g., [18]) studies the question of under 
which conditions fll.ip can be considered as a perturbation of the heat equation. The required 
regularity on b is usually expressed as b e B, with B a suitable function space. Typically the 
borderline spaces for which one can still prove most of the deeper resultsj are scale-invariant 
under the scaling (11.31) of b, that is, = ||6||e (see, e.g., [18~ | I29 T) @. The reason for this 

is as follows. The arguments in the proofs of the "deeper properties'!] typically have to work 
on all (small) scales and we therefore need to control b on all scales, which naturally leads to 
the scale-invariant spaces. 

Similar considerations can be made for elliptic equations of the form 

- Au + b ■ Vu = , (1.4) 
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with the elliptic scaling 

u(x) -> u (x \x) = u(Xx) , (1.5) 

b(x) — >■ b^ x \x) = Xb(Xx) . (1.6) 

Let us now consider the condition div 6 = and its consequences. (The relevant references 
include, for example, j36j [31] in the parabolic case and [22} [15] in the elliptic case.) Among 
the most important consequences are the following^. 

(i) The energy identity 

\u(x,t 2 )\ 2 dx + / / | Vu| 2 dx dt = / \u(x, ti)\ 2 dx (1.7) 

Jti JR n JR n 

is exactly the same as for the heat equation. 

(ii) The integral L„ u(x, t) dx is conserved: 

u(x,t 2 )dx= / u(x,tx)dx. (1.8) 

JM. n 

J. Nash showed in his famous paper [28] (inequality (8) on page 936) that one can obtain 
from (i) and (ii) the point-wise upper bound 

\G(x,t;y,s)\ < ■j-j—^y^ (1-9) 

on the fundamental solutions G(x,t;y, s), with C depending only on the dimension. There- 
fore, this bound also holds for solutions of (11. ip when div 6 = 0, with practically no other 
assumptions on b. The heuristic behind this estimate is that in an incompressible fluid, mix- 
ing can enhance the decay of, say, a temperature field but it cannot slow it down. Nash's 
simple argument proving this heuristics is very elegant. There are many other results in this 
direction, see for example 30. 3 1 . Bound (11. 9p can also be integrated in time to obtain 
(global) estimates of sup x \u(x)\ for the elliptic problem 

- Au + b-Vu = f , (1.10) 

with / G L n /2+s, a divergence- free b, and practically no other assumptions. 

Since the condition div 6 = has such strong consequences for the Loo-bounds, it is natural 
to ask about its effects on other properties of the solutions. For instance, can the standard 
assumptions on the drift term 6 needed, say, for the Harnack inequality be relaxed when 
div 6 = 0? Similar questions have been considered, for example, in [3TI [36] . 

It turns out that the condition div 6 = can be used to relax the regularity assumptions 
on 6 under which one can prove the Harnack inequality and other results. However, the 
effects are not as dramatic as in the case of Nash's upper bound (jl.9p . for which not much is 
needed beyond div 6 = 0. In particular, it seems that even with the condition div 6 = one 
cannot significantly "break the scaling" . Indeed, to be able to prove the "deeper regularity 
properties" of the solutions (as discussed above), we still need to assume that 6 belongs, 



To derive these consequences, one needs to assume that the formal integration by parts used to obtain 
them is valid. We are ignoring this technical issue for the moment. 
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at least locally, to a scale-invariant space B. The norm can be weaker than in the absence 
of the assumption div 5 = 0, but it still has to be scale invariant or stronger on the small 
scales. For example, the results of [31] imply that the Harnack inequality, Holder continuity 
of solutions, and the Aronson estimate for fundamental solutions^ remain true when b E 
LoofX^ 1 ), where L^ 1 denotes distributions which are first derivatives of bounded measurable 
functions. This should be compared to the condition b G L nj00 , which naturally comes up 
when the assumption div 5 = is droppecB Note that both L n and L^ 1 are scale-invariant. 

The assumption div 6 = can be used to reformulate equation (11 .ip in the following way. 
When div b = 0, we can write b = div d for an anti-symmetric tensor d = (<iij)0- Moreover, 
by introducing a suitable "gauge condition"!!, we can assume that the derivatives of d have 
similar regularity as b. Since b has dimension —1 with respect to the natural scaling of (11. ip . 
the tensor d has dimension 0, that is, it scales as 

d(x, t) ->■ d(Xx, \ 2 t) (1.11) 

when u is scaled by (II. 2p . 

Replacing b by the potential d, equation (II. ip becomes 

d t u - div(AVw) = (1.12) 

where A = I + d. This is a divergence-form equation with a non-symmetric leading term. 
Such equations (including the versions with lower-order terms) have been studied in [3T] 
under the assumptions that the coefficients are bounded measurable functions satisfying 
the ellipticity condition 

(A£)-£>H£| 2 - (1.13) 

The results of [31] show, roughly speaking, that most of the results which are valid for 
symmetric A are also true in the non-symmetric case. The transformation of (II. ip to (I1.12p 
has been used in many other works (see, for instance, [8]). 

In the elliptic case, Mazja and Verbitsky [22] studied (among other things) the bi-linear 
form 

(u,v) / (AVu) ■ Vv dx . (1.14) 



The form is obviously continuous in H l when A is bounded, but it turns out that the 
boundedness of the coefficients is not a necessary condition for the boundedness of the form. 
The form is still continuous on H l if the symmetric part of A is bounded and the anti- 
symmetric part of A is in the John-Nirenberg space BMO (bounded mean oscillation). This 
is a consequence of the following two facts: 

(i) If A is anti-symmetric, the form ( 11 . 141) can be factored through the determinants 3"'"") . 



5 ci(t - s)-" /2 exp[c 2 |z - y\ 2 /(t - a)} < G(x, t; y, a) < c 3 (t ~ s)- n ' 2 exp[c 4 |z - y\ 2 /(t - a)] , see ffl. 

6 Strictly speaking, as far as we are aware, when we do not assume div b — 0, most of the regularity results 
above are proved for b € L q ^ p with 2/p + n/q = 1 and p < oo (see [29]), but not in the borderline case 
p = oo, q = n. 

7 For n — 3, this corresponds to introducing the vector potential d such that b = curleL 
8 such as dkij + djkj + rf/j,fe = 0, which for n = 3 and b = curld corresponds to div d = 
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(ii) The determinants have "better than expected" regularity: when u,v are in H 1 , the 
determinants are not only in Li, but they are in fact in the Hardy space "H 1 , the dual space 
of BMO (see [5J). 

It is natural to expect that much of the classical regularity results for elliptic and parabolic 
equations with measurable coefficients in divergence form will remain valid if the leading part 
A is of the form A = a + d, with a symmetric, bounded and satisfying the usual ellipticity 
condition (I1.13p . and d anti-symmetric and belonging to BMO in the elliptic case, and to 
Loo(BMO) in the parabolic case. 

Indeed, let Q_ = lR n x E_ (with M_ =] — oo, 0[) and assume that 

A = a + d, (1.15) 

where a G L OQ (Q_; M nxn ) is a symmetric matrix satisfying 

ul < a < v- 1 ! (1.16) 

and d G /^(KL; BMO(M. n ; M™ xn )) is a skew symmetric matrix, that is, 

d = -d* (1.17) 

for all (x,t) G Q_. Here v > 0, I is the identity in the space M nxn of n x n-matrices and 
d* is the transpose of d. Let also B(x,r) be the ball of radius r centered at x G M n , and 
Q(z, r) = B(x, r)x]t — r 2 , t[ a parabolic ball in IR n+1 centered at point z = (x, t). Finally, let 
B = B(0, 1) and Q = Q(0, 1). We then prove the following parabolic Harnack inequality and 
Liouville theorem for suitable weak solutions (see Definition 12. ip to (11.121) . 

Theorem 1.1. If the matrix A satisfies conditions U.15\) - fil.l7\ ), then there exists C > 0, 
depending only on n, v, and ^d^L^i-x^BuoiS)), such that for any nonnegative suitable weak 
solution u to M.lty) on Q we have 

sup u(y,s)<C inf u(y,s), (1.18) 

(y,s)eQ(z R ,R/2) (y,s)&Q(z,R/2) 

whenever Q(z, R) C Q. Here, z R = (x,t — R 2 /2). 

Theorem 1.2. If the matrix A satisfies conditions U.15\) - fil.l7\ ), then the only bounded 
ancient suitable weak solutions to U.1S\) on Q_ are the constant functions. 

Remark. Of course, the corresponding elliptic results follow immediately by taking time- 
independent solutions. In addition, in Section 3 we provide a second — short and elementary 
- proof of the Liouville theorem for weak (sub) solutions (see Definition 13. ip to (11.41) in IR 2 . 

Recall that the norm in the space BMO{Vt; M nxn ) is 

BMO(n-M nxn ) = SU P I 31 J \d~ [d)x,r\ dx : B(x,r) <e Q 

^ B(x,r) 

with [d] X)7 . the average of d over B(x,r). 

We note that the space BMO is invariant under the scaling (11. lip , and hence these results 
are again in line with the argument that to preserve the "deeper properties" of the solutions, 
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one cannot "break the scaling". One of the goals of this paper is to present some evidence 
for this based on studying the failure of Liouville theorems under appropriate conditions. 

Let us first look at (11.41) in R n . By the Liouville theorem for (\1A\\ we mean the usual 
statement that a bounded solution in R n has to be constant. This is of course true for 6 = 0. 
For the time being let us assume that the vector field b is locally smooth, hence the solutions 
u are also locally smooth and the only obstacles to the validity of the Liouville theorem are 
global. 

The results of Stampacchia [34J imply the following: 
(L) IfbE L n (M. n ), then the Liouville theorem for A1.4]) holds. 

This is easy for n = 1, and for n = 2, there is also a relatively simple proof based on the 
energy estimate. The proof for n > 3 can be accomplished by using the Holder estimate or 
the Harnack inequality (see Sections 7 and 8 of [M]). If n > 3, then by Theorem 2.3 in |29j . 
the condition on b can be weakened to liminf^oo sup^^ 1 16| \L n (B(x,R5)) < c n for some 5 > 0, 
where c n > is a fixed dimension- dependent constant. This result implies in particular that 
(L) remains true when 

C 

\b(x)\ < -, — r for large \x\. (1-19) 
\x\ 

It is not clear to us whether condition (II. 19ft is sufficient in dimension n — 2, but this should 
not be a hard question. 

In dimension n — 1, condition (11.191) is sufficient when C < 1, as one can check by direct 
integration. With C > 1, however, f II . 1 9 j) is no longer sufficient. This can be seen from the 
example 

2x 
1 + x 2 

which was pointed out in this context to one of the authors in 1997 by Joel Spruck. The 
trivial extension of this example to higher dimensions is 

( 2xi \ 

b(xi, . . . , x n ) = I 2 j 0, . . . , ) and u{x\, . . . , x n ) = arctan(xi) . (1-21) 

We note that the vector field b in (11.211) belongs to the space (_BMO)~ 1 (M"), since 

ll ' ^log(l + x 2 ) (1.22) 



b(x) = — — - and u(x) = arctan(x) , (1-20) 



1 + x 2 dx 

and log(l + x\) G BMO(R n ). 

This example and Theorem 11.21 which establishes the Liouville theorem for b G (BMO) -1 
and div6 = 0, together show that the divergence-free condition can play an important role 
in Liouville theorems for equations with drift terms. On the other hand, we now provide a 
counter-example to the Liouville theorem with a divergence-free b on M 2 which is in some 
sense not too far from (BMO) -1 . Recall that the stream function of a divergence- free vector 
field b on R 2 is H : R 2 R such that 

b(x) = V^Hix) = (H X2 (x),-H Xl (x)). (1.23) 
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We therefore have 

- A + b ■ V = -div (A V), (1.24) 
where A(x) = I + d(x) has skew-symmetric part 

Theorem 1.3. There exists a divergence-free vector field b £ C°°(IR 2 ) with all derivatives 
bounded and a stream function satisfying \H(x)\ < Cm \x\ In In |x| for some C and all large 
enough \x\ such that 11 1.4ft has a non-constant bounded classical solution. 

This illustrates, to some degree, the important role of scale invariance of the assumptions 
in the Liouville result. In particular, it seems unlikely that one can significantly "break the 
scaling" even if we assume that div 6 = 0. 

We conjecture that similar negative conclusions can be arrived at when considering ques- 
tions about Holder continuity of solutions of f ll.4p (as well as the Harnack inequality). For 
example, it seems unlikely that the condition div b = is sufficient to get a C a -bound on 
solutions u in the unit ball B = 5(0,1) under the assumptions \u\ < C and ||&||z, n _ 4 < C. 
(Here we assume that all the functions involved are smooth, but only the indicated quantities 
are controlled, and we are interested in an a-priori bound.) 

Related to this are our last two main results, concerning distributional solutions u (see 
Proposition 14. ip of fll.4p in B with divergence-free b £ L\(B). The first establishes a loga- 
rithmic modulus of continuity of such solutions in two dimensions, depending only on ||6||li(b) 
and || w || /^(B). However, due to the low regularity assumed on the vector field b and u solving 
f ll.4p only in the distributional sense, our result is restricted to those solutions which can be 
obtained as weak-star Loo-limits of solutions with drifts in L,2(B). 

Theorem 1.4. Let B be the unit ball in M? and let (b m , u m ) £ ^(-B) x L 00(B) be a sequence 
of divergence- free drifts b m and distributional solutions u m to ljl-4\ ) with b = b m . Assume that 
u m are uniformly bounded in B and 

b m -> b in Li(B), 

u m ^u in L ^(B). 
Then the function u is a distributional solution to jl-4ty - Moreover, 

ueHUB)nC loc (B) 

and at the origin u has the modulus of continuity 

■ f , ^ ^(i + IHImb)) 172 ,, „ h 

sup \u(x) -u(0) < \\u\\ Loo {B) (1-25) 

iefl(o,r) V - log r 

with a universal C > 0. 

In three (and more) dimensions this result is false. Indeed, there exists no modulus of 
continuity of classical solutions depending only on ||6||li(_b) and ||w|| l x (b), and distributional 
solutions u £ Loo(B) fl ^(B) with divergence-free b £ L\(B) may be discontinuous. 



ON DIVERGENCE-FREE DRIFTS 



7 



Theorem 1.5. Let B be the unit ball in R 3 . 

(i) There is c > such that for each e > there is a smooth divergence- free drift b with 
\\b\\Li(B) — c and a smooth u with WuWl^ib) < 1? solving fll.4f) in B and satisfying 

u(0,0,e) -u(0,0,0) > c" 1 . 

(ii) There is a divergence- free drift b G Li(B) and a distributional solution u G H 1 (B) D 
Loo(B) of (11. 4p m i? which can be approximated by a smooth sequence (b m ,u m ) in the sense 



of Theorem \l-4\ but u is discontinuous at the origin. 



Our paper is organized as follows. In the next section, we develop local regularity theory 
for parabolic operators (11.121) under the assumption that the skew-symmetric part of A is 
in BMO, and prove Theorems 11.21 and 11.11 The important step of our approach is a higher 
integrability of suitable weak solutions. This allows us to adopt Moser's method for proving 
the Harnack inequality that implies Holder continuity of suitable weak solutions and Liouville 
type theorems for ancient suitable weak solutions. All these results hold true for the heat 
equation with a drift b G L^BMO' 1 ) as a particular case. In this connection, we would like 
to mention the recent paper [TO] , of which we learned while writing the present manuscript. 
In [TO], among other questions, the Cauchy problem for the heat operator with the drift term 
from L^BMO -1 ) has been considered and Holder continuity of solutions has been proved. 
The authors of [10J follow the Caffarelli-Vaseur approach [4J. In Section 3, an elementary 
proof of a Liouville theorem in the two-dimensional elliptic case is provided and Theorem 11.31 
is proved. Theorems 11.41 and 11.51 are proved in Section 4. 

Acknowledgement. GS was partially supported by the RFFI grant 08-01-00372-a. The 
other authors were supported in part by NSF grants DMS-1001629 (LS), DMS 0800908 (VS), 
and DMS-0901363 (AZ). LS and AZ also acknowledge partial support by Alfred P. Sloan 
Research Fellowships. 

2. Some results for parabolic equations 

The main goal of this section is to prove Theorems 11.21 and 11.11 We consider (11.121) in 
Q_ = R™ x R_, with the matrix A satisfying (jl . 15j) — fll . lTj) . We will study the so-called 
suitable weak solutions to (j!.12p . In what follows we will use the abbreviated notation 

B(r) = B(0,r), B = B(l), Q(r) = Q(0,r), Q = Q(l), 

as well as z = (x, t). 

Definition 2.1. Function u is said to be a suitable weak solution to equation (11.121) in the 
parabolic ball Q(R) if it satisfies 

ueL 2tOO (Q(R))nW^ (Q(R)), (2.1) 
J ud t tpdz= J {AVu)-V^dz Vcp eC^{Q{R)), (2.2) 

Q(R) Q(R) 
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and for a.e. t e] — R 2 , 0[, the local energy inequality 

to 

- J ip(x, t )\u(x, t )\ 2 dx + y y (pVu ■ aVudz < 



B(R) -R? B{R) 

U) to 



1 

) 

-R 2 B(R) -R? B(R) 



< 2 / / l M l 2 ^^ ~ I I (AVu) ■ Vyudz (2.3) 



holds for all non-negative test-functions p e C£°(B(R) x] — i? 2 , i? 2 [) . 

The function w : Q_ — > R is called an ancient suitable weak solution to HI .12)) . if it is a 
suitable weak solution to ( 11.121) in Q(R) for any R > 0. 

It is not clear whether one can show that any solution to (I1.12p . subject to assumptions 
(12. ip and (12. 2p . satisfies local energy inequality (12. 3p . In this respect the situation is similar to 
the Navier-Stokes equations: there is a certain cancelation due to the skew symmetric matrix 
d which works well in global setting, i.e., when initial-boundary value problems are under 
consideration. The corresponding procedure is relatively routine and leads to the existence 
of global solutions which satisfy the inequalities in Definition 12.11 at least locally. 

We now outline the main points of our approach. The structure of equation f 1 1 . 1 2 [) admits a 
modification of the technique developed by J. Moser in [25] [27] and get Holder continuity of 
suitable weak solutions. This property, together with scaling invariance, leads to the Liouville 
theorem. The main tool of proving Holder continuity is the Harnack inequality. We prove 
the Harnack inequality for smooth solution by the method of J. Moser. Extension of the 
Harnack inequality to suitable weak solutions is provided by higher integrability of the the 
spatial gradient. Here, our arguments use an approach due to M. Gianquinta and M. Struwe, 
see H2]. 



2.1. Local set-up and higher integrability. Equation (I1.12p is invariant with respect to 
translations and the following scaling 

u x (x, t) = u(Xx, X 2 t), A x (x, t) = A(Xx, X 2 t) (2.4) 

for any positive A. This allows us to reduce all considerations to some canonical domain, say, 
toQ = Q(l). 

So, we consider equation f 1 1 . 1 2 [) in the unit parabolic cylinder. Matrix A is split into 
two parts as in (TTTS]) with matrices a e L 00 (Q;M nxn ) and d G L 0O (— 1, 0; BMO(B; M nxn )) 
satisfying conditions (jl.ISp and (I1.17P 

In what follows, we shall denote by c positive constants depending only on n and v. We 

let NIUooC-BMO) = ||d|Uc«(-i,0;BMO(B)) and denote mean values by 

^ x ' r = Wr)\ j f ^ dy > ^ Z0 ' r = IQTrjT / u ( z ) dz - 

B(x,r) Q(zo,r) 

The main result of this subsection is the following theorem. 
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Theorem 2.2. Assume that u is a suitable weak solution to M.12\) in Q and matrices A, 
a, and h satisfy conditions U. 15\) - [1.17\ ). Then there exist two positive constants p > 2 
and C depending only on n, v , and ^^^(bmo) su °h that u G L P (Q(R)) for any R g]0, 1[. 
Moreover, the following estimate is valid: 

(tt^tt f \Vu\ p dzY < c( , * r. f \Vu\ 2 dzY (2.5) 
\\Q(R)\ J 1 1 / " \\Q(6R)\ J ' ' / K J 

Q(z ,R) Q(z ,6R) 

for all Q(z , QR) C Q with QR < dist(x , dB) and t - (QR) 2 > -1. 

This theorem is a consequence of the reverse Holder inequality, see [12] for further ref- 
erences. To prove the reverse Holder inequality, we need a Caccioppoli's type inequality. 
To formulate it, let us introduce additional notation. Fix a non-negative cut-off functions 
(p G Cq°(B(2)) and Xo{t) with the following properties: 

<p(x) = 1 x G B, x(t) = t < -4, 

Xo (t) = (t + 4)/3 -4<t<-l, Xo(t) = l t>-l. 
Now, for a point z = (xq, to) an d for R > such that Q(zo, 2R) G Q, we let 

Xt ,2R(t) = Xo((t - t )/R 2 ), <Px ,2r(x) = <p((x - Xq)/R). 

And then we can introduce a mean value of u as in [T2] 

u t0 ,2 R {t)= J u(x,t)(p 2 X0!2R (x)dx( J <pl 0t2R {x)dx^ . 

B(x ,2R) B(x ,2R) 

In our particular situation, we have 

Lemma 2.3. (Caccioppoli's type inequality) Under assumptions of Theorem \2.2i the follow- 
ing inequality is valid: 

to 

\u(x,t )\ 2 tp 2 X0)2R (x)dx + v J J X 4 2 ,2^o,2/?|Vu| 2 dz < 

B -IB 
to to 
<\j J \u\ 2( P 2 x ,2R 9 tXt ,2Rdz ~ J J Xt Q ,2R( aV ^) ■ ^vlo^Rudz (2.6) 

-IB -IB 
to 



J J Xt ,2R(( d ~ l d U,2R)^u) ■ Vcp 2 Xo>2R udz, 



-1 B 



where 



u(x,t) = u(x,t) - u X0)2R (t). 

Inequality \2. 6]) holds for a. a. to g] — 1, 0[, for all Xq G B, and for all R > subject to the 
additional condition Q(z ,R) C Q. 
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Proof. There are two important points to note. The first one is that for any skew-symmetric 
matrix oIq, depending on t only, we have 

d V« ■ Vipudz = (2.7) 

whenever ip £ C^°(Q). The proof is straightforward integration by part. 
The second point is that, see [12] . 



^ Ol MeLi(-l,0). (2.8) 
To see this, we take as test function in (12.21) the function (p 2 :Ql 2R( x ) r l('t) an d conclude 

d t u X(h 2R(t) = - J A{z)Vu{z)-V^ 2 X(h2B (x)olx/ J <pl ^ R (x)dx (2.9) 

B(x ,2R) B(x ,2R) 

Next, we replace u(x, t) with u(x, t) + u X0j2R (t) in local energy inequality (12. 3 1) and take 
ip = x 2( fx 2R with x from Cq(— 1, 1). Then terms which do not contain spatial derivatives 
can be transformed as follows 

\ J *o) + u XOj2R (to) \ 2 vl ,2 R {x)dx = 

B(x ,2R) 

= \ J \u{x,t G )\ 2 p 2 x ^ 2R (x)dx + ^\u X0) 2 R (t )\ 2 J <pl , 2R dx, 

B(x ,2R) B{x ,2R) 



and 



to 

1 

2 




( Px ,2R( x )\^( x i t ) + u x ,2 R (t) | d t x (t)dxdt 



-1 B 

1 

2 

-1 B 

to 




Vxn^iifc) l«<X *) l 2 9tx 2 (*)^ dt- 



X 2 {t)Ux Q ,2R{t)d t U X0: 2R{t)dt J (pl ,2Rdx+ 

B(x ,2R) 

+ ~|« a *,2Jl(to)|V(to) J ¥ 2 x ,2RdX- 



B(x ,2R) 

Now, the local energy inequality, together with the last two identities, implies 

to 

X 2 (t )\u(x,t )\ 2 dx + u J J x 2 (t)cpl 0!2R (x)\Vu(x,t)\ 2 dxdt < 

B(x ,2R) -1 B 
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to to 

<\J J vloi R (x)Hxit)\ 2 dtX 2 {t)dxdt- J J X 2 AVu ■ Vcp 2 Xo2R udx dt- 

-1 B -IB 
to 

2 



X (t)u X0 , 2R (t)d t u X0j2R (t)dt j <P Xo , 2R dx- 

-1 B(x ,2R) 
to 

-J J X 2 AVu-Vtp 2 X0t2R u X0! 2 R dxdt. 

-1 B 

By the (12. 9p . the sum of the last two terms is zero and from (12.71) it follows that 

i J \u(x,t )\ 2 x 2 (t )ip 2 Xo ^ R (x)dx+ 

B{x ,2R) 



to 

+V 



J j \ 2 {t)v 2 X0>2R {x)\Vu{x,t)\ 2 dxdt< 



-1 B 
to 



< \ I / rt , 2R (x)\u(x,t)\ 2 d tX 2 {t)dxdt- 




2 , 

-1 B 



to 

X 2 aVu ■ V<#£ 2R udx dt— 




-1 B 

to 



J J X 2 (d-[d} X0! 2R)Vu-Vif 2 X0j2R udxdt. 



-1 B 



Here, 



[d\x Mt) = \b(2R)\ / d< y x ^) dx - 

B(x ,2R) 

So, inequality (I2.6P follows if we choose the cut-off function \ in an appropriate way. □ 

Proof of Theorem \2. 6 A Using known simple arguments, we can derive from (12. 6p the following 
estimate 

to 

! = \J \u{x,t )\ 2 (p 2 XO!2R {x)dx + J J X 2 t ,2 R vl ,2 R \^u\ 2 dz< 

B -IB 

- C (^2 j \u\ 2 dz + ^ J (\Vu\cp xo>2R Xto,2 R )\u\\d - [d] X0j2R \dz 



Q(z ,2 R ) Q(zo,2H) 
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We now fix an arbitrary number s e]l, 2[. Let us denote as usual s' — s/(s— 1). Then the 
right hand side of the latter inequality can be estimated with the help of Holder inequality 
by 

to 

c r 2 , c 



— y |u| dz + — j ^ y |d- [dU^I* dx) x 

Q(z ,2ii) t -(2R) 2 B(x ,2R) 



X ( 

S(x ,2R) 



Applying Holder's inequality one more time, we find 

c 

Q(z ,2R) 

C 1 /" \ *7 

+— i?^ess sup sup ( i D , OD x| / |d- [c^^'dxj x 

-K t -{2R) 2 <t<t B(x ,2R)cB K \-D\ lJri )\ J J 

B(x ,2R) 

to 2 _ s 

>< y ( y iv«i Vxo,2iex t 2 ,2fl^) 2 ( y iwi^cte) 2s < 

t -(2R) 2 B(x ,2R) B(x ,2R) 

~F / ^ 2 ^ + ^^ll rf ll L -( BMO )( y l W lV' ,2ijXt ,2i?^) 2 X 
Q(z ,2R) Q(z ,2R) 

to 

r r 2— s 

x ( / f / \u\2^dx) s dt 2 



t -(2R.) 2 B(x ,2R) 

Summarizing our efforts, we have 

to 



i y |M(x,t |Vx ,2fl(a ; )^ + y y X^^/jl^ 2 ^^ 



-1 B 

to 



<c( S )(i+r 2 )R^ 2 y ( y \u\^d X y a dt, (2.10) 

t -(2R) 2 B{x ,2R) 

where T = ^^^(ba/o)- Now, let us discuss simple consequences of (I2.10p following [12J. By 
Poincare-Sobolev inequality, we have for 

s < (2.11) 

n — 1 



ON DIVERGENCE-FREE DRIFTS 13 

the following inequality 

( J <c(s)R n ^ +2 - n j \Vu\ 2 dx. 

B(x ,2R) B(x ,2R) 

Combining (12. lip and (I2.10p . we find 

J \u{x,t )\ 2 cp 2 X0i2R {x)dx < c(s)(l + T 2 ) J \Vu\ 2 dx. 

B Q(z ,2R) 

Hence, assuming that Q(z ,3R) C Q, we have the second estimate 

ess sup / |m| 2 (x, t)dx < c(s)(l + T 2 ) / \Vu\ 2 dx. (2.12) 

t -R 2 <t<t J J 

B(x ,R) Q(z ,3R) 

Now, our aim is going to be the so-called reverse Holder inequality. We first assume that 
the number s satisfies the condition 

2n 

Ks<- -, n = 2,3,... (2.13) 

2n — 1 

Obviously, fTXT3|) implies f l2ll|) and 

~" <1<^_< 2 , n = 2,3,... (2.14) 



2n-l _ 3~3n-2 

It is not difficult to show that under assumption ( 12 . 13[) there exist numbers < A < 1 
< n < 1, and 1 < r < 2 such that 

2s , nr 

2A + fi 



2 — s n — r 

X + H = 1 

nr 2 — s 

-\i = 1. 



n — r s 

Using these exponents, we derive from (12.101) 

to 



J \Vu\ 2 dz<c(s)(l + T 2 )R^-^ J ( J \u\ 2X+ ^dx) s °dt< 



Q(z ,R) t -{2R) 2 B(x ,2R) 

to 



< c ( s )(i + r 2 )R^-^ J ( J \u\ 2 dx) 7 \ j \u\^-r dx ) 7 "dt. 

t -(2R) 2 B(x ,2R) B(x ,2R) 

The last multiplier can be estimated with the help of Sobolev's inequality 

/ \Vu\ 2 dz < c(s)(l + T 2 )R { y~ 1)2 ess sup ( / \u(x, t)\ 2 dx) 2 x 

Q(z ,R) B(x ,2R) 
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r 1 
xR^( / \Vu\ r dz y 



Q(z (h 2R) 

To estimate the first multiplier on the right hand side of the last inequality, one can apply 
( I2.12p in the following way 



' dz 



\u(x,t) - u X0>2 R(t)\ dx < c J \u(x,t) - u Xo ,m(t)\ dx 

B(x ,2R) B(x ,2R) 

< c (s)(i + r 2 ) J \Vu\ 

Q(zofiR) 

for a. a. t &]t — (2R) 2 ,t [. Combining the latter inequality, we arrive at the reverse Holder 
inequality 

1 / |v^<c W( i + r V (^ / |v«|^)*x 



\Q(R)\ 

Q(z ,R) Q{z ,6R) 

' ; |Vu| r d^ 



\Q(2R)\ 

Q(z ,2R) 

which holds for some r g]1, 2[ and for any Q(z , 6R) C Q. This leads to a higher integrability, 

see ma. □ 



2.2. Moser iteration. To avoid some technical difficulties, we will assume that matrices a 
and b and solution u are sufficiently smooth in Q. Later we shall show how to remove this 
assumption. We also assume that our function u is strictly positive in the following sense 

u(z) >a R >0 \/zE Q(R) (2.15) 

for any < R < 1. Sometimes assumption (12.151) is not necessary, but for simplicity we will 
assume it is satisfied. We fix the following notation 

£ 2 (m) 

and, assuming that condition (12 . 13[) holds, let 



1 

2m 



2(n + 2) 

P = 



71 



2s p 

Q = 7> » 7 = - > 1. 

2 — s q 



Lemma 2.4. For any mi > m > 1/2 and for any < g < r with Q(z , r) C Q, we have 

sup u [z) < 

zeQ(z ,s) f r _ g\ 



sup u^{z)< Cl{n ^ S ^\ I »"""(:),/-. )-. (2.16) 



i 

Q(zo,r) 



where e = e(m ). 
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Proof. Set w = u m . For any m ^ 0, we can derive from fll . 12j) 

1 f 2m — 1 f 

- / ip 2 d t \w\ 2 dx H / ib 2 aS7w ■ Vwdx = 

2 J m J 

B(x ,r) B(x ,r) 

= -( y aVu; ■ wVtp 2 dx + J dVw ■ wVif> 2 dx) , (2.17) 

B(xo,r) B(xo,r) 

with a cut-off function ip satisfying: 

*fj(x,t) = <p(x)x(t), 
(p(x) = 1 xeB(x ,q), ip(x) = x B(x ,r), 
0< V <1, |V^|<^, 

*(t) = t<t -r 2 , x(*) = l *>*o-V, 

x(*) = - ~ r ° ~ 2 r ) ^o-r 2 <t<t -f? 2 - 

Next, we introduce the following sequence of exponents 

k = q, h = l%, 1 = 0,1,..., (2.18) 

If we let 

m i = l i m 1 /p, i = 1,2,..., 

then we have 

m i q = l i _ 1 m 1 , £ 2 (m;) = — 1 > £q, i = l,2, .... (2.19) 

Letting m = mi in (12.171) and taking into account (12.191) . we find 

sup / W *, t )|^ + * / *>|V«,|'*< 

U)-Q 2 <t<t J J 

B(x ,g) Q(zo,r) 

< ° J if) 2 \w\ 2 dz + cu- 1 J ip\Vip\w\Vw\dz (2.20) 

Q(zo,r) Q(zo,r) 



+c J \d - [d] X0)r \%l)\Vip\w\Vw\dz. 



Q(zo,r) 

The same arguments as in Section 2 show that the latter inequality gives us: 
\ w \l,Q{z , 8 ) = SU P / \w(x,t)\ 2 dx + I \Vw\ 2 dz< 

t -g 2 <t<t J J 

B(xo,p) Q(zo,q) 



(r - g} 



<^^(l + r 2 )r^( J \w\<dz 

Q(xo,r) 
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with s satisfying condition (I2.13p . By the known embedding theorem, see [18] . we have 
IMIj>,Q(*d,tf) - c \ w h,Q(zo,o) witn P = and > hence, 

(vim I i-H'sc( S , £o)( i + r)(-^-)x 



Q(zq,Q) 



71 1 f 1 

r\ a / 1 



It is worth noting that, under assumption (12. 13[) we have p > q. 
Our further steps are routine. We let 

R R 

= -Ri = o" + oi+T, r = Ri-i, i = l,2,..., 

in fl2~2T|) and find 

1 



\Q(Ri)\ 

Q(zo,Ri) 



u\ mih dzY l < 



Q(zo,Ri-i) 

for i = 1,2, .... After iterations, we arrive at (I2.17P with g = R/2 and r = R. General case 
is deduced from this particular one with help of known arguments. □ 

To see what happens if < m < 1/2, we have to introduce additional notation 
Q + (z ,R) = B(x ,R)x}t ,t + R 2 [, Q + {R) = Q + (0, R), 
Q(z ,R) = B(x Q ,R)x}t - R 2 ,t + R 2 [, Q{R)=Q{0,R). 
Lemma 2.5. For any < mi < 1/2 and for any < g < r provided Q(zo,r) C Q, we have 

sup ^)< C2(W iP ( / u m «{z)dz)\ (2.22) 
z&Q(z , e ) (r - g) <* v _ J J 

Q(zo,r) 

Proof. We replace the function \ with the following one 

X {t) = t>t' + r 2 , X(t) = l t<t' + g 2 , 

X(t) = ~ t+ I°V 2) t' + 2 <t<t' + r 2 , t' Q = t - f 2) V. 



r 2 _ Q 2 " - u V4 

Then from (I2.17p . we can derive (an analog of (12.201) ) 



sup / \w(x, t)\ 2 dx + e 2 {m)v I il) 2 \Vw\ 2 dz < 
:t<t +o 2 J J 



t <t<to+£? 

B(x ,g) Q+(4> r ) 
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<-, — °—, - / ip 2 \w\ 2 dz + as' 1 / ip\Vt/j\w\Vw\dz (2.23) 
{r - QY J J 

Q+(z' ,r) Q+(z' ,r) 

+c J \d- [d} XOjr \ip\V^\u\Vu\dz, z' = (x ,t' ). 

Q + (*'o> r ) 

Next, it is not so difficult to check that there exists a natural number k with the following 
property 



1 . , ,u..xA 
—mi < 

T 

And then, for this number k, we have 



mi < m'i = 7 ^ k < mi. 
7 2 2 



m[ < m 2 < ... < m' k < ~ < m' k+1 < 

where 

m: = ^ = 7 4 -x = y- fc -^, . = i,2,.... 

p 2 
and numbers U is defined by (I2.18p . It is easy to check that 

e 2 K)>7^-l, i = l,2, 
and then repeating derivation of (I2.2ip for w = with the same indices i, we find 



(|«bi / H^)W,r)( 



Vr — £ 



x 



0(4. e) 
-r\f 

V V|g(r) 

Q{*' ,r) 

Now, we consider (12. 24ft for 



X (-V(— 1^ / (2.24) 



and find 



r 1 r 

r — Ti, o — r»_i, n = — | : 

K * L ' 1 4 4 2* 

/ 

Q+(z' ,n) 

for % = 1,2, fc. After exactly iterations, we have 

V|Q(3r/4)| 7 1 1 / MQ(3r/4)| y 1 1 / " 

Q+(4,3r/4) Q(z ,3r/4) 
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1 



£c(s ' r) fe J lur " dz 

Since m' fc+1 > 1/2, we are in a position to apply Lemma 12741 letting there mi = m = m' fc+1 
and conclude that 

sup u<+i(z) < c(s,r)(— - / w^+^cfey. 

Q(2;o,3r/4) 



Taking into account definition (12.181) of and combining the latter inequalities, we find 

r i(i+7- fc )^/ i r 
sup v»*(z)< c(s,T) / 



and, by Holder inequality, we have 



sup w mi (^) < c(s,T)(— - / |ur'i 9 ^) 8 . (2.25) 



2GQ(z ,r/2) 



We may shift in time this estimate and show that 



i r — 
sup u m i(z)<c(s,T)(— / lup'^V, (2.26) 
zeQ+(z ,r/2) \r n + 2 J ) 

where 4' = (s ,*o) and *o = *o - From f l2T25|) and f[2?26]) . estimate f l2?22]) with e = r/2 
follows. General case is deduced from this particular one with help of known arguments. □ 

Lemma 2.6. For any e > and for any < g < r such that Q(z , r) C Q, we have 

sup u (z) < 

zeQ(z ,g) ( r _ £,) 



sup : r ' f,U/ -:: n ( / //-""I.:)*/:)". (2.27) 



9 

Proof. We let u = u~ £ and observe that by ( 11. 1211 . the function v satisfies 

d t v - divAVv < 0. 

We can repeat the proof of Lemma 12.41 with mi = 1 for v instead of u and then show 
2.3. Estimates oflnii. 

Lemma 2.7. Assume that u i a sufficiently smooth positive solution to equation U.lfy) and 
Q'(z ,R) = B(x , 2R) x]t — R 2 , t + R 2 [c Q. There exist two constants c 4 = c 4 (n, v) and a R 
such that 

ca R n + 2 

\{zeQ + (z ,R): -\nu-a R > s}\ < — , (2.28) 

s 

„ pra+2 

\{ze Q(z ,R) : -\nu-a R < -s}| < — . (2.29) 
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Proof. To simplify notation, we shift and scale our variables in the following way 

u R {y, s) = u{x + Ry, t + R 2 s), A R {y, s) = A{x + Ry, t + R 2 s) 

for (y,s) E Q' = 5(2) x] — 1,1[. Since equation (I1.12p is invariant with respect to this 
transformation, we may reduce our considerations to the cylinder Q' and, after proving our 
result for this particular case, get all the statements of the lemma with the help of inverse 
translation and dilatation. Without ambiguity, in what follows, we drop upper index R in 
the notation of functions u R and A R . 
So, if we let v = In it, then by (11.121) 

d t v - div(AVv) + Vv ■ aVv = (2.30) 

in Q'. Take and fix a smooth nonnegative cut-off function ip — ip{x) so that ip — 1 in B and 
ij) — outside -6(2). Multiplying equations (12.301) by ip 2 and integrating the product in x 
over -8(2) and in t over the interval ^2 [, we find 



vip 2 dx 



ta fa 
fa , , 




Vip z ■ AVvdx dt+ I I fV« • aVvdx dt = 

B(2) ti B{2) H -8(2) 




+ / / fVi) ■ aVvdx dt < 



and thus 

fa 

vijj 2 dx 

B(2) ti B(2) 

fa 

<cj J i)\Vip\\Vv\{\a\ + \d- [d] 0)2 \)dx dt. 

ti B(2) 

After application of the Cauchy-Schwartz inequality, we have the following estimate 

fa 

2 fa V 

vip dx 



<2 + 7T / / i> 2 \Vv\ 2 dxdt < c(r)(t 2 -*i). (2.31) 




fa 2 

B(2) fa B(2) 

From this point we essentially repeat arguments of J. Moser in [25], see Lemma 3 therein. 
We do this just for completeness. As it is pointed out in [25], we can choose our cut-off 
function ip so that the following Poincare-type inequality takes place 



\v{x,t) -V(t)\ 2 ip 2 {x)dx < c J \Vv{x,t)\ 2 ^j 2 (x)dx, 

B{2) B(2) 

where 

V{t)= J v(x,t)ip 2 (x)dx( J ip 2 (x)dx^j 1 . 

B{2) B{2) 
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Making use of this inequality, we can derive from (I2.3ip the following relation 
^(£ 2 ) - + 1 J J \v(x,t)-V(t)\ 2 dxdt<c 5 (n, v, r)(t 2 -*i) 

t\ B 

which can be reduced to the differential form 



^(0 + ^^ \v(x,t)-V(t)\ 2 dx<c 5 . 



r 

-(t) + cj 1 / \w(x, t) - W{t) \ 2 dx < (2.32) 



B 

One may make this inequality homogeneous with help of the shift 

w(x, t) = v(x, t) - 1/(0) - c 5 t, W(t) = V(t) - V(0) - est. 
This give us the inequality 

dW 
~~dt 

and the initial condition 

W(0) = 0. (2.33) 
For < t < 1 and s > 0, we introduce the family of sets 

B+(t) = {x G B : w(x, t)> s }. 

As it follows from ( 12 .321) and (I2.33p . for those values of parameters t and s, we have w(-,t) — 
W(t) >s- W(t) > on B+(t) and, hence, 

^(t) + c- 4 1 \Bt(t)\(s-W) 2 <0 

or 

The latter identity can be integrated and, as a result, we find 



1 



£4 
s 



which implies the first estimate ( 12. 28ft . Other estimate ( 12.29P can be established in the same 
way. □ 

2.4. Harnack inequality. 

Theorem 2.8. Let u be a positive sufficiently smooth solution to $1.12\) in Q. Then for any 
Q(zq,R) C Q we have the inequality 

sup u(z) < cg(ri, v, s, r) inf u(z), (2.34) 

z£Q(z 0R ,R/2) z€Q(z ,R/2) 

where z 0R = z — (0, R 2 /2). 
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Proof. Using translation and dilatation similar to those described in Section 4, one may 
consider our problem in a canonical domain, say, in Q'. Now our aim is to make use of 
estimates proved in Sections 3 and 4 plus some iteration technique in order to find a particular 
version of the Harnack inequality. It can be extended to the general case of Theorem 12.81 
with the help of covering methods. So, in this section, we follow [27] with minor changes. 
By Lemma [2 .7\ see (I2.29p . we know that 

\{zeQ: -In it -a < -s}\ < — (2.35) 

s 



for some constant a. As in [27], we introduce the following function 

<f(r) = sup In w(z) 

z€Q(zQ,r) 



with z = (0, 1/2), r e]l/2, 1/V2[, and w = e a u. 

Now, our aim is to show that there exist a constant 5 > 2 depending only on n, is, s, and 
■& E [1/2, l/y/2[ such that 

(p{&) < 5. (2.36) 
To this end, we derive from (12. 35ft the estimate 

w p dz < e p ^ r) -^r + e^ (r) (2.37) 
(p(r) 

Q(zo,r) 

being true for any positive mi — p/q and for any r >d. Let us choose m\ so that terms on 
the right-hand side of (I2.37P contribute to the sum equally. This suggests the following value 
for mi 

mi= — --ln-^. 2.38 
qp(r) 2c 4 1 ; 

Obviously, there exists a constant 5\ > 2 depending only on s and C4 such that if 

(p{&) > 5 U (2.39) 

then mi defined above belongs to the interval ]0, l/2[. If not, then 5 = 5±. If (12.391) holds, 
estimate (12.221) of Lemma 12.51 give us the relation 



mnp(g) < In ° 2 ^ + - In ( / w p dz) 

(r — o) 9 Q v J 7 



(r - g)' 

Q(zo,r) 

for any d < q < r. Recalling the choice of mi, we find from (I2.37P that 

w p dz < 2e^ {r) 

<9(zo,r) 

and thus 

1 c 1 

<p(0) < — In + ~v(r). 

I'll ^ r _ q Z 
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The latter inequality can be rewritten with the help of (I2.38j) in the following way 



<p(g) < 2^( r ) 



ln(c 2 (r - g) i ) 1 
miip(r) 



In c^r — g) 



-(n+2) 



l n Mr) 

111 2c 4 



Then one can consider two cases. In the first case, 



lnc|(r- g)- {n+ V < 1 
In Mil - 2 



2c 4 



and thus 



In the opposite case, we have 



/ii(n, i/,s,r) 



(r - £) 2 (™+ 2 ) ' 

Combining both cases, we find the following basic inequality 

Vte ^ \^ T) + (r - g)W + 2) 

for any 1/2 < ?9 < g < r < l/y/2. It can be iterated in the known way, see [27], and the 
result of these iterations can be expresses in the form 

So, f)2.36p is proved with 5 = max{#i, 5 2 }. 

Next, let z* = (0, 1) and v = u^ 1 . Then as it follows from Lemmata 12.61 and 12 .7\ see (I2.27P 
and QX2BD , 

c 3 



sup f (z) < 

zeQ(z ,£>) 



[r — g) i 



v q (z)dz 



Q{z ,r) 



and 



|{zeQ(^,l): \n(e~ a v) > s}\ < —. 



The same arguments as above show that there exists a constant 83(71, v, s, T) such that 

sup ln^iT 1 ^)) < 5 3 . 

This estimate, together with (12.36[) . implies a particular version of the Harnack inequality 

sup u(z) < e 5S3 inf u(z). 

The general case can be obtained from the particular case with the help of covering technique, 
see E3, Lemma 4, and translation and dilatation. □ 
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2.5. Nonsmooth case. 

Proof of Theorem Without loss of generality, we may assume that 

VuEL P0 {Q) (2.40) 

for some po > 2. To provide (I2.40p . we can apply Theorem 12.21 and scaling. Obviously, one 
can construct smooth approximations of matrices a and d with the following properties: 

a (jf) -»• a in L P (Q) 
d^ d in L P (Q) 

for any p > 1 and 

a a.e. in Q 
cr 3 ) — >■ d a.e. in Q. 
Moreover, matrices d^\ and ef^ satisfies conditions f ll . 1 5[) — ( fl . 1 Tf) with the 

same constants and \\d^\\ Loo (BMO) < \\d\\L x {BMO)- 

Then we consider the following initial boundary value problem 

d t w ij) - div(A (i) W j) ) = f ij) , 

W {3) \d'Q = 0, 

where d'Q is a parabolic boundary of Q and 

f {j) = <9 t «-div(A (i) Vu). 

Our claim is 

/^->0 in L^-l^F- 1 ), 

where if 1 is the completion of smooth compactly supported in B functions with respect to 
the norm ||w||2,b + ||Vu|| 2 £. Indeed, it is not difficult to show that 

ll/^lkc-x^) < (j{\a-a^\ 2 + \d-d^\ 2 )\Vu\ 2 dz)\ 

Q 

So, by (12.401) . the right hand of the latter inequality goes to zero. On the other hand, for 
wV\ we have global energy estimate 

Iw^Q^cWf^h^H^ 

which, in turn, means that 

\w®\ 2>Q ^0. (2.41) 

Now, we let = u — Obviously, is a unique solution to the following initial 

boundary value problem 

d t v U) - div(A {j) Vv u) ) = 0, 
(v U) -u)\ dlQ = Q. 
We know that possesses the following global properties 

fa® E L 2 (-l, 0; H- 1 ), W j) E L 2 {Q) 
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and, moreover, it is nonnegative on the parabolic boundary of Q and smooth inside Q where 
the equation for v^' can be reduced to the form 

- {a^vf\ k - d^vf = 0. 

Here, comma in lower indices stands for the differentiation with respect to the corresponding 
spatial variable and summation over repeated indices running from 1 to n is adopted. As 
it was shown in [18] , see Chapter 3, Theorem 7.2, therein, for functions satisfying equation 
above, the maximum principle holds and thus remains to be nonnegative everywhere 
inside Q. Obviously function v^' + j satisfies all the conditions of Theorem 12.81 and, hence, 

sup v^'(z) < c 6 inf v^\z). 

zeQ(z 0R ,R/2) Q(z Q ,R/2) 

Passing to the limit as j — » oo and taking into account f )2.4ip . we arrive at (11.18!) . ^ 

2.6. Liouville Theorem. In this subsection, we assume that u is an ancient suitable weak 
solution to equation (j!.12p which means that it is defined on Q_ = M n x] — oo,0[ and is a 
suitable weak solution in all parabolic balls Q(zq, 1) with zq = (xo,to) for any xq G R n and 
for any to < 0. Since our equation is invariant with respect to translation and usual parabolic 
dilatation, such a solutions will be suitable in all parabolic balls of the form Q(a) for any 
positive a. Now, we shall show its Holder continuity provided it is bounded. 

Lemma 2.9. Let u be an ancient suitable weak solution to equation M.12fy . Then there 
are two constants C7 and a which depend only on n, v, s, satisfying condition Iji2.13\) . and 
T = IMIUocfBMO) such that 

\u(z) - u(z )\ < c 7 \z - z \ a par sup \u{z)\ (2.42) 

zeQ- 

for any z and z from Q- with the parabolic distance \z — z \ par — \x — xo\ + \t — to| 5 - 
Proof. We let 

M R = sup u(z), M R/2 = sup u(z), 

z€Q(z ,R) z£Q(z ,R/2) 

rriR = inf u(z), rriR/2 = inf u(z). 

z£Q(z ,R) z&Q(z ,R/2) 

If we let v (z) = Mr — u(z), then will be a nonnegative suitable weak solution to equation 
(I1.12p in parabolic ball Q(z , R). By translation and parabolic dilatation, we can derive from 
Theorem 12.81 the following inequality for v 

inf v(z) = M R - M m > -{M R - u{z)) (2.43) 

z€Q(z ,R/2) C 6 

for all z G Q(zqr, R/2). On the other hand, for the same reason, we may apply Theorem 12.81 
to function w = u(z) — itir and find 

rriR/2 - m R > —{u(z) - m R ) (2.44) 
c 6 
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for all z G Q(z 0R , R/2). Adding (I2.43j) and (I2.44p . we arrive at the inequality 

osc(zo,R) — osc(z , -R/2) > — osc(2;o,-R), (2-45) 

c 6 

where osc(zo, R) = Mr — ttir, or 

osc(zo, R/2) < $osc(zo,R) 
with $ = 1 — 1/cq. After simple iterations, we have the series of the inequalities 

osc(>o, R/2 k ) < & k osc{z ,R) 
which can be reduced to the form 

osc(z , q) < c 7 g a, dosc(zQ, R). 
The latter is true for z G Q- and all < g < R < +oo and certainly implies (I2.42p . □ 

Proof of Theorem\TM We let M — sup 2g g_ |m(z)|. If we scaled our solution u and matrix A 
so that 

u R (y, s) = u{Ry, R 2 s), A R (y, s) = A(Ry, R 2 y), 
then as it is easy to see u R and A R satisfy the equation (11.121) in Q_ and 

v = M = M R = sup |^(^)|, r = T R = \\d R \\ Loo{BMO ). 

zeQ~ 



By Lemma [2. 9[ we have 

\u R (e) - u R (0)\ < c 7 \e\ a P arM R 
for any e = (y, s) G Q-. Making inverse scaling in the latter inequality, we find 



\u(z) -u{0)\ < c 7 \z\p ar —M 

for any z G Q- and for any R > 0. By arbitrariness of R, we show that u must be a 
constant. □ 

3. An elementary proof of an elliptic Liouville theorem in 2D and a 

COUNTEREXAMPLE 

In this section we explore the Liouville theorem for (jl.4p in two dimensions. Assuming 
that the divergence- free vector field b is in the space (BMO)' 1 , we provide an elementary, 
short, and self contained proof showing that bounded subsolutions (and supersolutions) are 
constant. Afterwards, we construct a counterexample to such a Liouville theorem for a 
divergence-free vector field whose stream function is bounded by In \x\ In In \x\ (^ BMO) for 
large |x|. This construction shows that the hypothesis b G (BMO)' 1 is quite sharp. 

If b is a smooth divergence- free vector field on M 2 , then it has a stream function H : M? — > R 
as in (I1.23P and we have (ll.24p . This relationship between b and A allows us to introduce 
the notion of a weak solution for very singular drifts. 
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Definition 3.1. Let b be a divergence-free drift from MO" 1 ^ 2 ), that is, H G BMO(R 2 ). 
We say that a function u G ^^(K 2 ) is a weak subsolution to f ll.4j) in M 2 , that is, a weak 
solution to 

-An + 6-Vw<0 (3.1) 
if for any nonnegative test function v G C^°(IR 2 ) we have 

y {AVu) ■ Vv dx < 0. (3.2) 

Weak supersolutions are defined by reversing both inequalities. 

We note that, as mentioned in the introduction, the bilinear form in (13 .2D extends contin- 
uously to compactly supported v G H 1 and hence in (13 .2D one can equivalently consider any 
nonnegative compactly supported v G if 1 (M 2 ). 

Theorem 3.2. Assume that b G 5MO _1 (l 2 ) is divergence-free, and let u be a weak subso- 
lution to [T4\ ) in M 2 . If u is bounded then u is a constant. 

Remark. If the drift b is not too irregular, for example, b G Z^ioc^ 2 ), then distributional 
solutions to (13. ip can be defined for u G L 2i io C (IR 2 )- In this case, bounded solutions are in 
Hl oc (E?) and satisfy H3.2D automatically, as we will show in the next section. 

The proof of Theorem 13.21 is an immediate consequence of the following lemma. 

Lemma 3.3. Let u be a bounded weak subsolution of 

- div(AVu) < in M 2 (3.3) 

where A(x) = a(x) + d(x) with a symmetric and d skew symmetric. Assume that there are 
A, A > such that for any x, £ G M 2 we have 

(a(x)0-ti>m 2 , (3-4) 
IMk»<A, (3.5) 
\\d\\ B MO<A. (3.6) 

Then u is constant. 

Remark. Note that we actually prove a Liouville theorem for bounded subsolutions. This 
is only possible in two dimensions. In higher dimensions one needs a to be a solution in order 
to show that it is constant even in the case of the Laplace equation. 

Proof. Without loss of generality we can assume that u is nonnegative (otherwise we can add 
a constant). Let r] be the test function 



if \x\ < 1 

losR 



if Id > R 
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We take v = urj 2 in (I3.2p to obtain 

> J {AVu) ■ V(urf ) dx 

> J rj 2 (aVu) -Vudx + 2 J wq(AVu) ■ Vt] dx 

> A J | Vu\ 2 7] 2 dx + 2 J uri(AVu) ■ Vr] dx 

(all integrals are over IR 2 unless otherwise indicated). Therefore we have 
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(3.7) 



A / [Vu|V dx < 2 



urj(AVu) ■ Vr) dx 



(3. 



We need to estimate the second term. For the symmetric part of A, we have 



ur](aVu) ■ V77 dx 



< A||uVr/|| L2 ||r/Vu|| L2 

<l\\vvu\\l + c\\v v \\l 2 

< -llnVull? + — — 
- 8 M ' ML2 logR 



(3.9) 



for a constant C depending only on IHIl^, A, and A. 

Let k be the average of k in B(R), the disk of radius R centered at the origin. It is easy 
to check that 



urj(dWu) ■ \7i] dx = 



Now we estimate the contribution to the variable skew-symmetric part of the coefficients 
using Holder inequality: 



UT)(dVu) ■ V77 dx 



< C 



\ 



1/4 



/ \d-dA (|^w) 1/4 (/^i 2 ) 1/2 

\B(R) J 



< j I t/ 2 |Vw| 2 dx + C 



( 


\ 






\B(fl) 


J 



1/2 



« 4 |V7?| 4 



1/2 



(3.10) 



Since d is a BMO function, we have 



\d-d\ 4 dx < CR 2 



B(R) 
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and by direct computation using that u is bounded, 



J u 4 \Vr]\ 4 dx < 



C 



R 2 (\ogR) 4 



Therefore, ( 13 . lUj) gives 

w](dVu) ■ Vr] dx 



< t [v 2 \Vu\ 2 dx + — ? (3.1i; 
4 7 (log/?) 2 

Adding (13. 9ft and (13. lip , we estimate the right hand side of (13. 8ft to obtain 

r C C 

|Vti| 2 dx < \Vu\ V dx < 



hgR (\ogR) 2 

B{1) R2 

for a constant C independent of R. We conclude the proof by taking R — > 00. □ 

Proof of Theorem\Tj^ Let /i : Mq -> Mq be such that h(s) = e 1 ~ e s for s G [0,e e ] and 
h(s) = In sin Ins for s > e e . For x = (xi,x 2 ) define x = min{|xi|, \x2\} and 

H{x) = C sgn(x 1 x 2 )h(x) (3.12) 

with C large. If now H = rj * H for some radially symmetric smooth mollifier rj supported 
on the unit disc and b = V ± H, then the hypotheses of the theorem are satisfied. Moreover, 
if = {x J ± X2 > \x\ \ + 2}, then for some c > independent of C we have 

bi{x) = and sgn(x 2 )6 2 (x) > cCh'(x) if x G K + U K~ . (3.13) 

Let B t be the 2-dimensional Brownian motion with B = 0, defined on some probability 
space (f2, J 7 , P), and for x G M 2 let X t x = X^(u) be the stochastic process with Xq = x and 
satisfying the SDE 

dXf = b(X*)di + V2dB t . 
It is then well known (see, e.g., |30j) that if v solves 

d t v + b ■ Vv - Av = 

on R 2 with t>(0,x) = t>o(x), then 

v(t,x)=E(v (X?)). (3.14) 

We let v = 2x{x\x 2 >o}- Clearly v(t,x) G [0,2] by ( 13.141) . and (&i(xi, -x 2 ), b 2 (xi, -x 2 )) = 
(b 1 (x 1 ,x 2 ),-b 2 (x 1 ,x 2 )) gives 

v{t,x u 0) = 1 for all (t,Xi). (3.15) 

So d t v (t, xi, 0) = and obviously sgn(x 2 )d t v (0, Xi, x 2 ) < for x 2 7^ because of v(t, x) G [0, 2] 
and the choice of Vq. This and the maximum principle for d t v give sgn(<9 t v(t, x)) = — sgn(x 2 ) 
and so there exists w(x) = lim^oo v(t, x). Parabolic regularity shows that u is a (bounded) 
solution of (Oj) . 

Let = : (Xf(w)) 2 ^ for all a G [0,t]} and A x = [j t>0 A t>x . Then (EH, (EUD, 
and the strong Markov property for Xf imply 

v{t,x) = (1 + sgn(x 2 ))P(A,x-) + 1 - P(A,x), 
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u(x) = (1 + sga(x 2 ))¥(A x ) + 1 - ¥{A X ). 

We will now show that P(^4( Q)4 )) = P(A( _4)) > (the equality holds by symmetry), which 
implies u(0,4) > 1 > u(0, -4). 

The law of iterated logarithm (see, e.g., [30]) implies that F(A) > for 

A = {u : \B t \ < (3tlnlnt) 1/2 for all t > e e and \B t \ < 1 for all t G [0,e e ]}. 

For any u G A, let t > be the first time such that Xj ' 4 ^(u;) G dK + . Then clearly t > e e 
and so 

|(^i°' 4) )i| = |V^(S S ) X | < (Gtlnlnt) 1 / 2 (3.16) 
for all s G [0,*]. Since X s (0 ' 4) G X+ for these s, we have X s (0 ' 4) = Kxi 0,4 ^! and (l3~13l) gives 

for s G [0,t]. This means 

1 /2 

(^ ( °' 4) ) 2 > cC t + V2(B t ) 2 + A> (^-^(tlnlnt)^. 

If we choose C > 18c _1 , then this and (I3.16P give 

(Xi°'% > 2(6tlnlnt) 1 / 2 > (Gtlnlnt) 1 / 2 + (6e e ) 1/2 > KXi ' 4 ))^ + 9, 

contradicting X^ '^ G dK + . Therefore X^ ' 4 ^^) G K + for all t > and w6 A. This means 
that A C A (0 ,4) and so < P(^4 (0i 4)) = P(A( ,_4)). Hence u(0,4) > 1 > w(0, -4) and the 
result follows. □ 

4. On a modulus of continuity in 2D and a counterexample in 3D 

In this section we prove that distributional solutions of (11.41) with a divergence-free b in a 
two dimensional domain ftcK 2 are continuous with a logarithmic modulus of continuity that 
we estimate explicitly. The modulus of continuity depends on a local bound for the H l norm 
of u and it does not essentially depend on any quantity associated with the vector field b. If 
b G Lij oc , then for suitable u (see Theorem 11.41) we can estimate the modulus of continuity 
in terms of the norm of u instead of H 1 thanks to a local energy inequality. Because of 
the low regularity assumed for the vector field b, the a priori estimates are hard to extend to 
distributional solutions and this presents some technical difficulties that are explained below. 
The estimate is a version of the classical result that functions in the border-line Sobolev 
spaces which satisfy the maximum principle^ are continuous, with logarithmic modulus of 
continuity. 

We also show in this section that the same type of regularity result does not hold in three 
dimensions. Indeed, we construct an example of a function u G L^B) D H 1 (B) (recall that 
B = B(0, 1) is the unit ball) and a vector field b G Li(B), such that u solves ( II. 4p in the 
distributional sense and is discontinuous at the origin. 

9 Sometimes the terminology "monotone in the sense of Lebesgue" is used in this context, see e. g. [21] . 
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Proposition 4.1. Assume that the drift b £ L 2i io C (^) is divergence- free and u £ L tX3 (f2) 
solves fll.4p in the distributional sense (that is, 



J b ■ Vvdx = = J (uAv + bu ■ Wv)dx 



(4.1) 



for all v £ C^°(Q)). Then u £ Hi oc (Q) and if B(xo,r) d Q, there is a constant C depending 
only on r such that 

1 /2 

1 1 Vw| 1x2(5(2,0^/2)) < C (l + | \Li(B(xa,r))) 1 1«| 1^(5(^,0) 

Proof. The claim w £ i^ c (fi) is obvious from bu £ L 2 ,ioc(^)- 
This means that 

J b ■ Vvdx = = J (Vtt - 6u) • Vvcte (4.2) 
n n 

for v £ if 1 (f2) compactly supported in Q. Let 77 be a smooth bump function such that 

77 = 1 in B(xo,r/2) 

rj = in fi \ B(xo, r). 
We take t> = m/ 2 in ( 14. 2 p to obtain 

= J \Vu\ 2 r] 2 + 2wyVu ■ Vr/ — b ■ Vuurj 2 — 2b ■ Vr/u 2 ^ dx 

B(xo,r) 

I Vu\ 2 r] 2 + 2urjVu ■ V77 — 6 • Vr]u 2 rj dx, 



B(so,r) 

where we have used 

0= J b ■ \7(u 2 r] 2 )dx = 2 J b ■ Vr/u 2 !] + b ■ Vuurfdx 

B(x ,r) B(x ,r) 

with u 2 r) 2 £ H l (Q). Therefore 

I Vti| 2 r/ 2 dx = J — 2nr/Vn ■ Vr? + b ■ Vr/u 2 !] dx 

B(xo,r) B(x ,r) 



~\ J l Vu lV d x + J 2 ^ 2 |V^| 2 + b ■ Vr]u 2 7] dx 



B(xo,r) B(xo,r) 

and thus 

1 



2 

B(x ,r/2) B(x ,r) 



Vu\ 2 dx < j 2u 2 \Vr]\ 2 + b ■ Vr]u 2 r] dx 

B(xo,r) 

<C(l + ||6|| Xl (B( a ,o.r)))||«| 



|2 
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□ 

It is worth noting that all statements of Proposition 14.11 hold true in higher dimensions. 

We now find the modulus of continuity for functions satisfying the maximum principle and 
a bound in H 1 . Note that in two dimensions, the space H l is borderline with respect to 
the Sobolev embeddings to spaces of continuous functions. The monotonicity of oscqb(t)U is 
the extra assumption used in the theorem below to actually obtain an explicit modulus of 
continuity. 

Theorem 4.2. Let u G H 1 ^) and assume that for any r G (0, 1) the maximum principle 
holds in B{r): 

max-u = maxii, 

B(r) dB(r) 

min w = min u. 

B(r) dB(r) 

Then u satisfies the following modulus of continuity estimate at the origin 



C 

sup \u{x) — ~u(0)| < \\^ u \ 

x&B(r) V-logr 



L 2 (B) 



for some constant C independent of u. 

Proof. Let r G (0, 1). We want to estimate oscb^w = maxe( r ) u — min^) u. 

i 

|Vd 2 dx = / |Vd 2 da ds 



B\B(r) r dB(s) 

since |Vw| 2 = u 2 a + u 2 v where u a is the tangential derivative and u v is the normal one, 

i 

> / / \u a \ 2 da ds 



r dB(s) 

Rewriting the integral using polar coordinates (s6 = a), 

i 

= S~s I \ U0 ( sd ^ 2d6 ds - 

r dB 

Since H l (dB) C C a (dB) from the one dimensional Sobolev imbedding, 

i 

/C 
— (oscaB( s) w) 2 ds. 



32 ON DIVERGENCE-FREE DRIFTS 

From the maximum principle, oscqb( s )U is monotone in s, therefore, 

1 

/C 
— (osc d B( r )u) 2 ds = (-C logr)(osc dB (r)u) 

r 

Taking square roots of both sides we obtain 

OSC B ( r )M = OSC 9B r r \U < \\VU \\L 2 (B) 

V— logr 



□ 



Consider now a drift b G L\(B) and let u G L^B) be a distributional solution to ( 11. 4ft . 
We are interested in whether u is still a continuous function and, if so, how to estimate its 
modulus of continuity. We do not know the answer to this question. However, it is in the 
affirmative if u is an appropriate limit of solutions with L 2 drifts, as in Theorem 11.41 

Proof of Theorem \l-4\ The first claim is immediate from the definition of distributional so- 
lutions. Moreover, Proposition 14.11 and Theorem 14.21 show that u m are locally uniformly 
bounded in H 1 as well as locally uniformly continuous with the modulus of continuity from 
( ll.25p . and the second claim follows. □ 

Finally, we show that that Theorem 11.41 does not hold in higher dimensions in general. 

Proof of Theorem \1.5l As in the previous section, we will again consider vector fields with 
b(Rx) = Rb(x), where R(x\, x 2 , £3) = (xi,x 2 , —xs). In addition, b will be axisymmetric with 
respect to the X3-axis and with no angular component. Such divergence-free vector fields can 
again be obtained from a "stream function" H : M.^ x M. — > M. with H (0, z) = as 

1 
2p~ 2 



b H (x) = V x 



— (xiH^p, z), x 2 H z (p, z), -pH p (p, z)) , (4.3) 



2p- 

where p = \J x\ + x\ and z = x 3 . Notice that again we have bn • Vi? = = 6^- (x 2 , —X\, 0), 
so H is constant on the streamlines of bn, and bu has no angular component. 
We now pick a G (|, 1) and for p 2 + z 2 < 1 and p > we let 

{p 2 z- 2 p a < \z\, 

z 2 p- 2 \z\ a < p, (4.4) 

(pl^l)^ 1 -^ 1 ^ \z\ < p a and p < \z\ a . 

Finally, for some large C we define H = CH and b = b Ho . 

Notice that H is continuous and vanishes on the axes, and bo(Rx) = Rbo(x). We also have 

- b (x) = C-^-, for {x\ + x 2 2 ) a/2 < \x 3 \ (4.5) 



F3 



as well as b G Li(B) (because l&o^)! < c\x\ 2 for some c > due to a G (|, 1)). Moreover, 
H is smooth except on 

P = {p 2 + z 2 < 1 : p > and \z\ G {p Q , p 1/o , 0} } , 
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so b is smooth except on 



S=<xeB: [ Jxj + xl,x 3 ) e P> . 



We therefore let H £ be smooth such that H e (p, —z) = —H £ (p,z) and H £ = H outside the 
e-neighborhood of P, the vector field b £ = bn e is also smooth and |6 e (x)| < c|x| -2 on B, as 
well as 

gm||6 e -6o|k ( B) = 0. (4.6) 

Clearly V ■ b £ = for e > 0, so (14. 6[) gives V ■ b = in the distributional sense. 

For each e > we now construct (smooth) u £ using b £ in a way similar to our construction 
of u using b in the proof of Theorem 11.31 We are here on B so we set boundary conditions 
u £ (x) = sgn(2 3 ) on OB and thus consider the stochastic process 

dX^ £ = -b £ (Xt' £ )dt + V2dB t . 

with Xq' £ = x and stopping time 

r £ = M{t>0\Xt' £ edB). 

We therefore obtain 

- Au £ + b £ -Vu £ = (4.7) 

where 

u £ (x) = E (sgn ((X£*) 3 )) = sgn(x 3 )P ((Xf' e ) 3 ^ for all t G [0, r e ]) , (4.8) 

Each u £ is a smooth solution of (14 .7p and they are uniformly Holder continuous away from 
the origin due to |6 e (x)| < c|x|~ 2 . Therefore there is a sequence Sk — > and m such that 
u £k — > uq locally uniformly on B \ {0}. (In fact, 

\imu £ (x) = u (x) = sgn(x 3 )P ((^ t x ' ,0 ) 3 + for all t G [0, r ]) 
for any x G B, provided we set Uq(0) = 0.) But this, Hi^Hl^ < 1, and ( 14. 6 ft show that 

/(M Af + 6 m • Vf )<ix = lim / (-u £fc At> + b £k u £k ■ Vv)dx = 

S B 

for any t» G C oo ( J B). Thus 

- Am + b ■ Vm = (4.9) 

in the distributional sense. The proof of Proposition 14.11 applies to each u £k , implying that 
their weak limit u G Hl oc (B). Since for e > we have u £ (Rx) = —u £ (x) and ||6 e ||li(b) is 
uniformly bounded, we only need to show 

limlimuJ0,0,z) > (4.10) 

24,0 e— >0 

to conclude the proof of both (i) and (ii) . 
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In fact, let us consider instead of (0, 0, z) with z > any y G B with y 3 > and \fy\ + y% < 
|?/3^ a . Let K y be the cut-off cone 



(2ys) 



l/a 



ir y = jx G M J : ^/x 2 + x 2 , < x 3 j C R 2 x 

with upper and lower base consisting of discs D y , E y centered at (0,0, 2y 3 ), (0,0, \y 3 ) and 
with radii \(2y 3 ) 1 t a \ l |(2y 3 ) 1 / Q . Notice that its tip would be at the origin, were it not cut off. 
Then \Jx\ + x\ < x 3 on K y because a > §, so (14.51) holds on K y fl B. Let a be the exit 

time of Xt' £ from K y fl -B, which is the same for all e < y\ because b £ = bo on K y fl -B in 
that case. We will show that 

F(X y ' E G D y UdB)>l- e' y ^ +2/ \ (4.11) 

provided C from the definition of H is large. This is sufficient since X y,£ G D y U 95 means 
either 

(Xf £ ) 3 ^ for all t G [0, r £ ] (4.12) 

or (Xl<% = 2y 3 . In the latter case we have ((X y a > £ )\ + (X^) 2 ) 1 / 2 < \(X' y > £ fJ a , so we can 
bootstrap (14.111) and obtain (I4.12p with probability at least 

L- lo S2 s/3j o 

(1 - e -^ 3+Va ) > n (1 - e -^ 3+2/a ) ee m > 0. 

fc=l j=—oo 

Thus (14.81) gives M e (y) > m for any y & B such that y 3 > and \jy\-\-y\ < \y^ a and any 

£ ^ y 3 1/a - The claim of (OOP now follows immediately. 

It remains to prove (14. lip for some large C independent of y. The point here is that Xf' E 
starts well inside K y (specifically, dist(y, dK y ) > dyl^ a for some d G (0, 1)) and the (strong) 
drift —h quickly pushes it towards D y while the Brownian term will not affect this picture 
much during the short time needed to reach D y , at least with probability close to 1. We have 

X y / = y- f b Q {Xr)dt + B a , (4.13) 
Jo 

as well as 

P(|£ t | < dyl /a for all t G [0, 8C7- 1 ?/ 3 ]) > 1 - e~ y ^ +2/a , (4.14) 
provided C is large enough. Since the vector — bo(X y,£ ) points 'inside' the mantle of K y for 
t G [0, a) (because of (14 .5p and the fact that the cut-off tip of K y is the origin) and has 

a positive third component, and dist(y, dK y ) > dy^ a , this means that with probability at 

least 1 — e~ y 'i , the process Xf' £ cannot exit K y through the mantle or the bottom E y 
before time 8C" 1 ?/ 3 . But we have (-6 (Xf £ )) 3 > C{2y 3 )- 2 for t G [0,a], so (14331) . (Oil . 

and y 3 + C^y^C^yi - cy X 3 ja > 2y 3 yield 

-, „ -3+2/a 

F(X y > £ eD y UdB and a < 8CT 1 ?/ 3 ) > 1 - e~ v * 
This proves (14. 1 1 j) and the result follows. □ 
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